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If gamma-ray bursts (GRBs) produce high energy cosmic rays, neutrinos are expected to be 
generated in GRBs due to photo-pion productions. However we stress that the same process also 
generates electromagnetic (EM) emission induced by the production of secondary electrons and 
photons, and that the EM emission is expected to be correlated to the neutrino flux. Using the 
Fermi observational results on gamma-ray flux from GRBs, the GRB neutrino emission is limited 
to be below ~ 20GeV m~^ per GRB event on average, which is independent of the unknown GRB 
proton luminosity. This neutrino limit suggests that the full IceCube needs stacking more than 370 
GRBs in order to detect one GRB muon neutrino. The Fermi observations of GRBs also imply that 
the ratio between energy in the accelerated protons and electrons is fp < 10. 

PACS numbers: 95.85.Ry, 14.60.Pq, 98.70.Rz, 98.70.Sa 



/. Introduction. The sources of high energy, > IPeV, 
cosmic rays (HECRs) are expected to produce high en- 
ergy neutrinos by pion production. The detection of 
high energy neutrinos wiU help to identify the origin of 
HECRs. Gamma-ray bursts (GRBs) have long been pro- 
posed to be one of the strong candidates of extragalactic 
HECRs[ll-Ia|, and was expected to produce high energy 
neutrinos [3, [y]. Currently the IceCube, operating in the 
full scale, is the most sensitive TeV-scale neutrino tele- 
scope, and is believed to reach the level of GRB neutrino 
flux. The recent reported flux limits from IceCube in its 
uncompleted configuration have put stringent constraints 
on GRB neutrino emission The observations by Ice- 
Cube in full scale will soon give even more stringent re- 
sults. 

The comparison of the latest non-detection of GRB 
neutrinos by IceCube and the positive GRB neu- 
trino prediction challenges GRBs being the source of 
HECRsQfs^. However, the predicted flux depends on 
some uncertainties in the GRB model. First, the neu- 
trino flux is proportional to the proton luminosity from 
the GRB jet Lp, which is unknown in GRB model and 
regarded as a free parameter. Second, the neutrino flux 
is proportional to the fraction of proton energy that is 
converted into pions /tt, which is further dependent on 
other uncertain parameters, e.g., the bulk Lorentz fac- 
tor of the jet, r, and the radius of the GRB emission 
regions, i?em@, ^M- Given the fiducial values of these 
parameters one can calculate the GRB neutrino flux, as 
done by Guetta et al. [13] and IceCube0-[ll, but the flux 
is subjected to the uncertainties. For example, if the ra- 
tio between energy in accelerated protons and electrons 
fp = Lp/Le is taken to be fp = 10, the predicted neu- 
trino flux by P, [l^ challenge the current non-detection 
of neutrinos by the IceCube, but using fp — 1 even the 
photosphere model of GRBs, which has small Rem (hence 
large /^), survive the current detection limit [l6l|. A sys- 
tematic consideration of the parameter values makes rel- 
atively more reasonable prediction [l^l, but the result still 
suffers the assumptions of parameters. 



Here we investigate the neutrino production in GRBs, 
and emphasize that the neutrino flux could be related 
to the electromagnetic (EM) one in the pion production 
processes. Then we derive a constraint on GRB neutrino 
flux based on the Fermi observations of GRB EM emis- 
sion. As shown below, this normalization of GRB neu- 
trino flux to gamma-ray emission suffers much less model 
uncertainties, in particular, independent of the proton lu- 
minosity Lp. Without special mention, we will use in the 
following the convention Qx — Q/IQ^ and cgs units. 

//. EM and neutrino correlation. The neutrino emis- 
sion from GRBs are caused by photopion production. If 
HECRs are generated in GRB outflows, they can inter- 
act with the observed intense photon field and produce 
charged pions via p + j n + . The charged pions de- 
cay via the primary mode, 7r+ /i^ + f^, and muons fur- 
ther decay via /i+ e'^ + + So the initial neutrino 
flavor ratio in the source is : : =1:2:0, 
but due to neutrino oscillation the flavor ratio observed 
in distance becomes : : Ri 1 : 1 : ifUj. 
We can see that final products induced by one charged 
pion decay are four leptons, each roughly shares the pion 
energy, 

« ei.^ « « e^^ « K/4. (1) 

The generated secondary high energy electrons may eas- 
ily convert their energy into gamma-rays by radiation 
processes. Thus there should be a straightforward rela- 
tion between the muon neutrino flux and the gamma-ray 
flux induced by the secondary electrons, 

E.,.^Ff. (2) 

This means that the neutrino flux may be normalized to 
and then constrained by observed gamma-ray flux, which, 
as shown later, appears mainly in GeV scale. 

The secondary electrons are very energetic and cool 
rapidly by synchrotron radiation. Let us calculate the 
synchrotron photon energy. For a flat proton distribu- 
tion with index p w 2 {Epdrip/dEp cx Ep~P), and a 
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typical GRB photon spectrum with Band-function pa- 
rameters, a-y = 1 and (3j = 2, the energy distribution of 
the generated secondary pious is dominated by pions at 
^ E-^^i, — 0.2Ep h, where Ep^i, is the energy of protons 
that interact with photons with spectral-break energy e^, 
at Delta resonance [Sj], Ep^b = 1-3 x lO^^Fj ge^^^y^V, 
where eb^MoV = et/lMeV. Given that each pion-decay 
generated electron carries one fourth of the pion energy, 
the synchrotron photon spectrum emitted by the sec- 
ondary electrons is mainly dominated by photons above 

esyn ^ 6(eB/ee)^/'i52'At"^e,-LvTeV, (3) 

where L is the GRB gamma-ray luminosity. At is the 
variability timescale in GRB light curve, and and es 
are the fractions of jet kinetic energy that are carried 
by accelerated electrons and magnetic field, respectively. 
High energy photons may interact with the GRB MeV 
photons and produce electron/positron pairs, 77 — > e^. 
The 77 optical depth exceeds unity for photon energy 
above [illiSi 

e-y^ ~ max(lA^_2^l5L^2^ 0.3r2.5)GeV. (4) 

Thus Esyn 3> e-yy, i.e., the synchrotron photons emitted 
by secondary electrons may be easily trapped and con- 
verted into electron/positron pairs. This initiates an EM 
cascade: the pairs emit synchrotron photons; the pho- 
tons are converted back into pairs; the pairs emit pho- 
tons again. The photons can escape once their energy 
decays into below e^^ during the EM cascade. There- 
fore, the total energy of the pion-decay induced elec- 
trons finally show up as a EM cascade emission piled 
up around e^^, typically e^^ ^ O.lGeV— a few lOO's GeV 
for r ~ 100 — 1000. This is exactly the energy windows 
of Fermi/LAT. Therefore, we may expect to constrain 
GRB neutrino flux by using gamma-ray flux observed by 
Fermi/LAT. 

A. Additional gamma-ray components. It should be 
noted that the observed gamma-ray flux, F^, may be 
contributed not only by the secondary electrons induced 
EM cascade emission. Let us write F^^ ~ fF-y with 

F-f = F:^^ + F^'^'^ and / < 1. The additional contri- 
butions to the gamma-ray flux could be several origins. 
First, there may be leptonic component (F^), emitted 
by GRB accelerated electrons, other than secondary elec- 
tron/positron pairs from hadronic processes. Second, the 
produced neutral pions in p7 interactions also introduce 
additional gamma-ray component (F^ ), which will be 
considered in more details in the following. 

There is comparable possibility that the fry interac- 
tions produce neutral pions, p -f 7 ^ p + vr". Each neu- 
tral pion decays into two photons, tt*^ — > 77, and each 
photon carries one half of the primary pion energy. Ap- 
proximately we assume the neutral pions follow the flat 
energy distribution of the charged pions, thus the energy 
distribution of the secondary photons from neutral pion 
decays is dominated by those above 

e^,b « 0.5E^,b = 1.3 X IQi^F^ se.-j^^yeV. (5) 



This is well above the 77 absorption energy (eq|4]), so 
the pion-decay induced photons should be absorbed and 
lead to an EM cascade emission. However, for very high 
energy photons the Klein-Nishina effect will reduce the 
77 interactions. For photons with energy larger than^iol. 

eesc = 10i6Ai:ir-2L526ZLveV (6) 

the optical depth reduces to below unity and the photons 
can escape. Here eL = lei.kcvkeV is the low energy 
cutoff of the GRB spectrum, due to, for example, the 
synchrotron self absorption, which is important typically 
below keV scale [2l[. Thus the photons with energy of 
^Ti.b < e < Ceps are trapped inside the GRB emission 
region and their energy is transferred into EM cascade 
emission. 

On the other hand, the generation of neutrinos and 
secondary electrons at high energy ranges suffers from 
the fast cooling of charged pions and muons. Electrons 
are generated by muon decays, which are slower than 
their synchrotron cooling for muon energy higher than 

i?^,, = 3 X lOi6(ee/eB)i/'ri5Ai_2L5-2'^'eV0, M M- 
The secondary electrons that dominate the EM cascade 
emission are mainly produced by charged pions with en- 
ergy of E^^b < < (4/3)^^,^. 

For a flat primary proton distribution, the gener- 
ated secondary pions are also flat in energy distribu- 
tion, approximately E'^duTr/dET^ (x E'^. Assuming in 
the P7 processes, comparable energies are carried by 
the generated charged and neutral pions, the ratio be- 
tween the energy of secondary electrons from tt"*" de- 
cays to that of secondary photons from tt^ decays is 
~ ilog(|i^^,c/-E^,6)/log(ecsc/e7r,6), with the factor 1/4 
due to that electrons carry one fourth the energy of pri- 
mary pions. If the gamma-ray flux is only contributed 
by EM cascade emission from secondary electrons and 
photons, and given the relation of eq (O, the / factor is 
then calculated as 

Ff + Ff ^ \ hgilE^jE^.b) + log(eesc/e..fc) ' 

(7) 

Plugging the values of i?^,c, £'7r,b, Ccsc and e^^.b we have 
/ « 0.36 for typical parameter values. If there are other 
contributions to the gamma-ray flux, e.g., the leptonic 
component, it should be that / < 0.36. 

///. Fermi Limit to neutrino flux. The Fermi Gamma- 
ray Space Telescope provided a brand new window to 
GRB observations, especially for the GeV scale. We show 
here how low the GRB neutrino signal has been con- 
strained to be by Fermi observations of GRB gamma-ray 
emission. 

Fermi satellite has two instruments. Gamma-ray Burst 
Monitor (GBM), with wide field of view (FOV) and be- 
ing sensitive to MeV scale emission, and Large Area Tele- 
scope (LAT), narrower FOV and sensitive to GeV emis- 
sion, especially from 0.1 to few lOO's GeV. The LAT has 
detected roughly 8% of the GBM-triggered GRBs that 



3 



occurred within the LAT FOV. For these LAT-bright 
GRBs, the analysis by Zheng et al.jl^ has shown (their 
Table 2) the photon fluence (i.t., time integrated pho- 
ton flux) of each LAT detected GRB in the 0.1-300GeV 
range. We then get the av erag e photon flucence for one 
GRB, (/)bright = 108phm^^[3i|. Since there is quite few 
photons detected above few lO's GeV, the photon num- 
ber in 0.1-lOGeV range is practically equal to that in the 
0.1-300GeV range. In their recent paper[23| the Fermi 
team has analyzed the LAT-dark GRBs in roughly 3 
years operation, which are the other 92% GBM-triggered 
GRBs that are in the LAT FOV. The upper limits to the 
gamma-ray flux in 0.1-lOGeV range for each GRB has 
been listed in their Table 1. We use the upper limits in 
the last column to calculate the average upper limit of 
all these LAT-dark events, which is (/)dark = 34.6phm~^ 
(0.1-lOGeV range). Let us regard ^bright as the upper 
limits to the bright GRBs, then we can calculate the av- 
erage upper limit to all GRBs, including both bright and 
dark GRBs, as 

0iimit = 8% X ^bright + 92% X 0dark = 40phm-2 (8) 

at 0.1-lOGeV energy range. Given the photon fluence, 
one still need the spectral profile in order to know 
the energy fluence. Assuming a flat photon spectrum, 
e^dn-y/de oc e^~^ with 7 = 2, the average limit to the 
fluence of one GRB is 

F^,ii„,it = 20GeVm-2 (9) 

in 0.1-lOGeV range. For softer photon spectrum with 
7 = 3, the limit becomes i^^jimit = 8.5GeVm~^, smaller 
by a factor of 2.3, whereas for harder spectrum with 
7 = 1.5, F^, limit = 75GeVm~^, larger by 3.7. Given the 
correlation between neutrino flux and gamma-ray flux, 
the average neutrino fluence per GRB is F^,^ < /f 7, limit- 

IV. Neutrino detection rate. Besides the normalization 
of neutrino flux, one still need the neutrino spectral form 
in order to calculate the neutrino detection rate by ex- 
periments, like the IceCube. For a GRB spectrum with 
Band- function parameters and 13^, and for a flat pro- 
ton distribution with p = 2, the GRB neutrino spectrum 
generated by the p"f interactions can be approximated as 
dn^/de oc e^"" at e < ei, duv/de oc e~^'' at ei < e < £2, 
and dn^/de oc e at £2 < £ (For simplicity hereafter we 
take £ = e^^) where the spectral indice are = 3 — /J^, 
= 3 — a-y, and = j3i, + 2, and the break energies 
are £1 = E^^^b/A and £2 = E^^c/S. The normalization of 
the neutrino flux is obtained by, given F^,^ — fF^, the 
requirement that e{dni,/de)de = /i^-y, limit- For typ- 
ical GRB spectrum with = 1 and /3j = 2, we get 
the specific photon flux at £ = £1 is dn^/ds{e = £1) = 
/-F7,iimit/ei[(3/2) 4- log(£2/£i)]. Notc, the broken power 
law form is a good approximation to the spectral profile 
of GRB muon neutrinos resulted from a full numercial 
calculation with the effect of neutrino mixing, e.g. [25l|. 

The effect area of the IceCube in its 40 strings con- 
figuration and averaged over incoming angles can be ap- 
proximated accurately by a function: log(Aj.§^'^/cm^) = 



a + bix + 622;^ + b^x^ + h^x'^ where x = log(£/lGeV), 
a = -17.27, bi = 9.16, 62 = -1.67, 63 = 0.142 and 
64 — —0.00466. The effect area of the full scale IceCube 
(86 strings) is about 3 times larger, A^^^^ « 3^J^g4°[2i|. 

Now the average neutrino number that can be de- 
tected in one GRB cann be calculated as A'dot = 
Acs{dnij/de)de. For the neutrino spectrum with 
F^^ < 20/GeVm-2, £1 = 7 x lO^^eV and £2 = 
lO^^eV, the detection rate by the 40-string configuration 
is AriC40 ^ 2.5 x 10'^ f (for 7 = 2) in one GRB, so more 
than l/Nl^^° = 400// GRBs are needed to stack in or- 
der to detect one muon neutrino from GRBs. This GRB 
number are larger than that analyzed by the IceCube 
collaboration in ref 0\ . More GRBs are needed to be ob- 
served for positive detection. For the full scale IceCube 
with larger effect area, and for / — 0.36, one needs to 
stack at least 370 GRBs to detect one GRB muon neu- 
trino (or > 100 and > 850 GRBs for 7 = 1.5 and 3, 
respectively) . 

V. Conclusion and discussion. We show that in the 
P7 processes there is a correlation between the generated 
EM (electrons and photons) and neutrino fiuxes. The EM 
emission cascades down to typically GeV scale, before 
escaping from the GRB emission region, thus one can 
use the Fermi/LAT observations to constrain the GRB 
neutrino flux, which implies that for the IceCube in its 40 
string configuration, one needs > 400// GRBs to detect 
one GRB muon neutrino, or > 370 GRBs for the full 
IceCube (taking typical value / — 0.36). 

Since the gamma-ray flux can be contributed by lep- 
tonic processes, e.g., generated by radiation from GRB 
accelerated electrons, the LAT observations only give an 
upper limit to the EM emission from interactions, and 
then an upper limit to the neutrino flux. It should be no- 
ticed that it is possible that / ^ 0.36. 

The derivation of neutrino flux using the EM-neutrino 
correlation is less dependent on the uncertain parameters 
in GRB models. For example, unlike the previous calcu- 
lations of GRB neutrino flux, the calculation here does 
not need to assume the (non-thermal) proton luminosity 
in GRBs, Lp, or the ratio fp. Note that the derivation is 
not parameter free; the location where EM cascade emis- 
sion piles up depends on the bulk Lorentz factor F (eq. 
2]). However, observational results usually suggest that 

r - 102 - lo^fHllllll. 

GBM detects roughly ~ 10^ GRBs each year. The av- 
erage limit to the neutrino flux per GRB then suggests 
that the full IceCube detects < 3 muon neutrinos asso- 
ciated with GRBs each year. This is smaller by order 
of magnitude than the prediction by Q , who normalizes 
the neutrino flux to the observed flux of ultrahigh en- 
ergy, > lO-'^^eV, cosmic rays (UHECRs). This may chal- 
lenge the assumption that GRBs are the origin of the 
observed UHECRs. Otherwise, the local, within ^ 100 
Mpc, HECR production rate may be larger by one order 
of magnitude than in the distant, earlier Universe, since 
UHECRs come from within the GZK sphere whereas the 
neutrinos can reach the Earth from the edge of the Uni- 
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verse. 

Given the neutrino flux, one can constrain the HECR 
production in GRBs, i.e., the ratio between energy in 
accelerated protons and electrons. It is derived that the 
protons of energy Ep^f, typically loses a fraction /^^b ~ 0.2 
of its energy by pj interactions in GRBs [5], and that the 
ratio between luminosity in neutrinos and protons isfioj 



Lp 8 ' log(i?p_inax/ Ep^^in) 



1 , 1 

8 •' ' 5 



(10) 



Here i?max (-Emin) is the maximum (mini- 
mum) energy of accelerated protons. Thus 
Lp ^ 200Ly^. The proton to electron ratio is given 
by fp — Lp/L^,yicY = {Lp/ L^^){L^^/L^, Mcv) ~ 
200(F,^/i^^,Mov) < 40(F^,Mev/10-5ergcm-2)-i. In 
fact, the wide energy range observations by Fermi may 
have make a clearer constraint. The Fermi observations 
show that the ratio between flux in LAT and GBM 
is typically i'y,Gcv/-^7,McV ^ 0.1[29]. Noting that 
/ ~ L^^/L^^GeY < 0.36, by using fp ~ Lp/L^, mcV = 
(Lp/L^^)(L;,^/i^_Gov)(i7,Gcv/i7,Mcv), wc have 



fp<7 



f 



^7,GcV 



200L^„ 0.36 O.IL, 



7,McV 



(11) 



If the GeV gamma-ray flux is dominated by leptonic com- 
ponent, then / ^ 0.36. This indicates that the Fermi 
limit starts to challenge the use of the fiducial number 
/p = 10 as in the previous literatures. 

There may be some caveats that the Fermi Limit can 
be avoided in some cases. First, we have used the GeV 
scale, 0.1-few lOO's GeV, flux to constrain the neutrino 
flux, but it could be that the cutoff photon energy is 
much larger, e.g., e-y^ 100 GeV, so the observations in 
< 100 GeV range do not make sense. Second, if for un- 
known reasons it happened that the generated secondary 
electrons, i.e., generated by the charged pion decay or 77 
pair production, do not radiate at all, the EM-neutrino 
correlation would not exist. 

The calculations here assume that the p7 interactions 
happen in a region of radius R « T^cAt. There are 
other GRB models where this assumption may not hold, 
then one should remind that the neutrino flux is model 
dependent [11]. 
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